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1 Introduction 



The conventional point of view [1] connects quark confinement with a hnear potential which 
is expected for two static quarks inserted in the vacuum of pure gluodynamics. It is usually 
assumed that the colour field of the quarks interacting with non-perturbative gluon vacuum 
fluctuations shrinks into the tube. This leads to a linear growth of the energy with the 
distance between quarks and confines them. The situation could drastically change, how- 
ever, when light quarks are introduced in the theory. Here the string between static quarks 
instantly tears and created quark pairs screen the linear potential of the pure gluodynamics. 
This means that the specific properties of QCD with light quarks are determined by con- 
figurations of gluon fields rather different from the configurations providing confinement in 
pure gluodynamics. 

The simple example of the phenomenon discussed is the well-known Schwinger model 
presenting the quantum electrodynamics of massless fermions (quarks) in two dimensions 
[2,3]. Here the original Coulomb potential is indeed proportional to the distance but the 
interaction with massless quarks screens it out. As the result, one can observe in QED2 

i) the absence of charged (coloured) states [2-4]; 

ii) confinement of quark quantum numbers (chirality, fermion number [5], flavour for mul- 

tifiavour QED2 [4,6]); 

ill) the chiral vacuum condensate [3,6,7]; 

iv) the independence of quark jets and "soft" hadronization [8,9]. 

The above features of QED2 seem to be similar to that expected in QCD. 

In this paper we just investigate the electromagnetic (e.m.) fields which provide the above 
features of QED2. We show that these are fields for which the Dirac operator determinant 
D{A) vanishes. We divide the above fields in two types to be, respectively, A and B. 

The fields of the first (A) type are necessary generated in QED2 by charged systems. 
They are responsible for charge confinement in QED2. In order to prove this we consider 
the evolution operator S{T) of the system in the quark representation. This operator is the 
sum over all physical states 

^(T) = $:e-^^"^|vl>„><vl>„|, (1) 

h 

En being the energy of the state. Matrix elements of the operator S{T), Eq.(l), can be 
represented as functional integrals over the fields with boundary values fixed at initial and 
final momenta of time. In the Coulomb gauge which we shall use in the paper we have:Q 



S{T) = jDAoexp -i j ^(fx j D'^D'^expiJ "^{id + gAo)'$d^x 

I J ^,,,(0) '- 



(2) 



^This gauge in convenient since here are no unphysical degrees of freedom. This gives possibihty to 
present a clear physical interpretation of all phenomena considered 
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For QED2 we have constructed S(T), in ref. [10]. 

Matrix elements of the evolution operator determine the probability amplitudes to have 
any quark configuration in any physical state. The charge confinement means that such 
amplitudes must vanish for the configurations where some charge (quark q) is situated far 
from other compensating charges (antiquark q). We shall demonstrate that this is the case in 
QED2 and the amplitude to find quark on distance r from antiquark falls as r~^. Investigating 
the structure of the integral, Eq.(2), one can discover the reason for such a dependence: the 
self-consistent field A of the distantly situated quark annihilation the quark determinant 
D{A). This is the field of the class A. Since the amplitude is proportional to D{A) it 
naturally vanishes at large r. The explicit proportionality to D{A) appears in Eq.(2), after 
integration over fermion fields. We have [10]: 



Here Sext{Ao,t) ~ D{Aq) is the evolution operator in the given external field Aq (explicit 
form see in Sec. 2). Its matrix elements are determined by chosen quark configurations and 
are expressed in terms of the fermion Green functions in the external e.m. field Aq (Eq.(5)). 

Now we see that if for some quark configuration a self-consistent field Ao nullifies D{A) 
then this configuration is prohibited; it cannot appear in evolution of the system. It is 
this mechanism which forbids existence of charged states in QED2. As we have stated the 
field of every charged quark creates the field , which makes determinant to vanish. The field 
strength of the discussed field behaves as 1/x for large distances due to the effect of the 
vacuum polarization. 

The fields with D{A) —>■ always create new quark pairs. These pairs screen origi- 
nal charges. If the total screening is not achieved, a new configuration also corresponds to 
D{A) — > and the process repeats. It is well known that D{A) determines vacuum polar- 
ization effects. Hence, the described processes represent the mechanism of charge screening 
by means of the vacuum polarization. 

As it had been noted by many authors [6,11,12], the screening mechanism of QED2 
reminds in some aspects the well-known Higgs one. However, there is an essential difference 
between them. The QED2 models contain no charged condensate because any configuration 
with separated charges cannot be present in the evolution operator S{T). Therefore, there 
are no charged complexes in the physical vacuum. We had seen this fact investigating the 
quark structure of QED2 vacuum states [7]. Charge screening is not statistical (in average) 
as it happens in the Higgs case, but is the exact one. The wave functions of excitations 
are here the eigenfunctions of the charge operator with vanishing eigenvalues. It is just the 
spectrum of the type expected for the real world [1]. 

The fields of the second class (B) govern chiral properties of the model. They also make 
D{A) 0, but decrease with a space distance faster than x~^. They can be characterized 

^In fact, |i?(A)p is the probability of the process where no quark pairs are produced. 




(3) 
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by a non-zero value of the integral 

T T 

Qtop{T) --£11 d'xE{x,t) --£^1 J d^xe,.F^^{x), (4) 



which is the two-dimensional winding number [13]. The S-fields also inevitably produce 
gg-pairs. As we shall see 2|Qt„p| particles of these pairs necessarily have very small momenta 
p ~ being considered in the finite space volume V {V ^ oo). They are delocalized over 
the whole volume. Such particles do not contribute to any local quantities and are, in fact, 
unobservable. In the V ^ oo limit, the non-conservation of chirality K (which follows from 
the Adler-Bell-Jackiw (ABJ)-anomaly [14]) can be interpreted in QED2 just as a generation 
of delocalized particles with small momenta p ~ 1/V by B-ficlds. 

The 5-ficlds with integer Qtop play an important role in the formation of quark conden- 
sates. For Qtop ~ N exactly pairs qnqL (or qnqL for < 0) are generated by such a 
S-type field. The chiral condensates of QED2 with one and several fiavours [2,6] are formed 
by similar chiral complexes [7] (Sec.3). 

We shall see in Sec.3 that i?-fields prevent a change of chirality for developing quark 
configurations. Owing to this property the exchange of chirality between hadrons and the 
chiral condensate is impossible. As a result, chirahty of hadrons is conserved and appears 
to be the definite fixed quantity {K — Qi). 

In some aspects the properties of B-fields arc similar to those of instanton fields [15] in 
QCD. However there are essential differences between them, as it is discussed in Sec. 2. 

The second main topic of this article is the hadronization process investigated in Sec. 4. 
We consider there the production of a g^gi-pair by an external source. This, widely discussed 
in QED2, process [3,8,9] is an analogue of the e+e~*hadrons reaction. It is also governed by 
a field with D{A) — > (B-type). Such a field produces just that one gR^L-pair with finite 
momenta {p ~ m) which screens charges and chiralities of the both jets. Simultaneously, 
the same field creates another q^qL-^paSi with the vanishing total momenta which joins the 
vacuum condensate. 

Thus the screening of charges and chiralities for both jets is finished at the finite time 
t ~ {m — /"K is the hadron mass), when QtopiT) becomes of the order of unity. At 
larger times t > l/mwe have only neutral independent jets with K = 0. From these times, 
there begins a rapid growth (proportional to t^) in the number of quarks with momenta 
p ^ m. At times t ~ p/m'^ their distribution is transformed into the hadron distribution 
of the parton model (1/p). The leading pair of quarks loses its energy, and after a time 
t ~ Pin/n^'^ participates in the formation of the hadrons. The produced hadrons with p < pm 
become spatially separated from the leading quarks of the initial pair. Quarks exist outside 
hadrons only on the light cone. 

2 Properties of fields where D{A) 

EM fields which make the determinant D(A) correspond to those quark configurations 
which are prohibited for evolution and production. In order to prove this statement let us 
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consider the evolution operator S{T) for a finite time interval T. As it was shown in Ref.[10], 
the operator S{T) can be represented as a functional integral over A. In the Coulomb gauge 
(^1 = 0) we have Eq.(2) for S{T). 

The integration in Eq.(2) goes over all EM fields Aq in the time interval < t < T 
including end points t — and t — T. The evolution operator for massless fermion in an 
external field Aq is equal to [10] 

Sea^tiT) = D{Ao)ew\ [ dxdx' [4{^)GP{xT,x'0)ai{x') 

{•' i=R,L 

+ bi{x)GPix, 0; x', T)bl{x') - a+(x)G'f ^(x, T - e; x'T)ht{x') 

- hi{x)Gf\x,s;x',Q)ai{x')\}; (5) 

£ ^ +0. 

Here D[Aq) is the determinant of the Dirac equation in the external field Aq\ o,rl{^)j ^rl(^) 
(a~ = a) are the creation and annihilation operator for right (R) and left (L) quarks and 
antiquarks: 

^i{x) = ai{x) + bt{x), i^R,L, = ^(1± 75)^(0;). (6) 

The Green functions Gi in an external field have the form [10] 

gP{x, t; x', t') = Gf\x, t; x', t') exp J ig j dti j dxiAQ{x, t) 

I 

X Gf\x,t-x',M)-Gf\x\t'-x^M)\], 

Gf\x,t]x',t') = [27ri{t - t' T {x - x) - iO sign (t - t'))]-\ (7) 

g is the dimensional coupling constant of QED2. 

For any processes, matrix elements of S{T) are determined by quark configurations in 
the initial (aj(x), and final {af{x),bf{x)) states. These matrices appear in Eq.(3) as 

coefficients in front of corresponding numbers of a^, 6^-operators. To obtain matrix elements 
we must expand the exponential in Eq.(5). Then the integral, Eq.(3), will be of the Gauss- 
type one. Its value can be obtained substituting the saddle-point field Aq, in every term of 
the integrand. 

Hence, if D{Ao) — > 0, the configuration discussed does not contribute in the norm of any 
physical state vector: 

*(r) = 5(r)*o- (8) 

It does mean that this configuration cannot appear as a result of evolution of the system. If 
charged states exist, then distantly situated charges would be represented in matrix elements 
of the 5'(T)-operator with a finite probability. There are no charged states when all such 
matrix elements decrease with space distance r between charges. As we shall see, in QED2 
all the above fenomena are due to fields providing D{A) — > 0, at r ^ 00. 
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In this section we discuss general properties of the fields D{A) — > 0. We take here fields 
as external ones. Self-consistent fields of QED2 considered in the next section, are of the 
D{A) type, too. 

The Dirac determinant for massless fermions in an external field has the following form 
[10]: 

Ini^(^) --^J^j dtdt'E{p,t)E{-p,t)^^. (9) 

' ' 

— g'^/fT, E{x,t) = —dAo/dx is the field strength in the Coulomb gauge. 
All the field configurations discussed correspond to finite field strength. Then D{A) can 
vanish only because of the divergence of the integral in Eq.(9) at p ^ 0. Therefore, we shall 
consider two types of the field D{A) — > 

(A) E{p,t) ^ signp/i(i), (B) E{p,t) ^ Mt), p^O, (10) 

where fi{t) are some functions of t. These fields have difi'erent asymptotics for x 00. The 
field strength for the A-type fields decreases as 1/x, while E for the B-type ones decreases 
more rapidly. 

Evaluating a divergent part of the integral, Eq.(9), we obtain the following formulae for 
the D{A) using space volume y as an infrared regularization: 

D{A) = 1/VP'^^\ PiT) Tdt [ dxV.P. I MM (11) 

27r 7o J J X — y 

for the A-type fields and 

D(^) = l/yl«-^(^)l, Q,,,{T) = -^£dtJdxE{x,t) (12) 

for the B-type. Qtop is the "topological charge", of the field (see Sec.l). 

The main contribution to the norm of state vectors for the fields D{A) comes from the 
particles with small momentum p ~ 1 /V. Indeed, since the norm of wave functions in Eq.(8) 
is conserved, the small factor D{A) should be compensated by the contribution of quark 
configurations expanded over the large volume V. Particles with small momenta p ~ 1/V 
represent just such configurations. So the probability of no pair creation in D{A) fields 
is suppressed by the factor ~ 1/y to be compared with the probability of configurations 
with quark pairs. 

In order to make explicit the role of quarks with small momenta for A and B fields, let 
us consider the number nR{p,t) of quarks produced by such external field in the fermion 
vacuum: ^o{t) — exp(—iHext{t))\0 >. The simplest way to calculate nii{p,t) is to use the 
bosonization procedure. For the fields of the B-type these calculations have been performed 
in detail in Sec.3 of Ref. [9]. To find ^o(^)) one must solve the Schrodinger equation in an 
arbitrary Aq external field taken in the Coulomb gauge. For this aim we express ^'(a;), ^'"'"(a;) 
operators by means of bosonization forms [9] and calculate nii{p,t) as follows: 

n«(p,t) = (*o(^)l4(p)«ii(p)l*o(^)) = Jdxdye'P^^-y^ 



5 



dxdy e'P'-'^-y^ 
X exp{27ri[Q;(x, t) — a{y, t)]}. 



(13) 



We use the normalization condition to be {a'^{p), a,{p')} = V6p^pr. In this case we have 

dk e-*'^^*-^) 



27r k 



a{x,t) — j dti j dxiQ{t — X — ti + xi)E{xi,ti) = j 



^R{k,t) = g /271 J dhj dxie'''^^'-'''^E{xi,ti). 
For small k we have for A and B types respectively: 

$R(/c,t) = isignkpit), = -Qtop{t). 

Substituting Eqs. (14) and (15) to Eq.(13) and taking into account the relation 



My 

we obtain the following formulae: 



1 _ r dk 

-x-iO) ~ y 27r' 



-ik{y—x) 



(14) 



(15) 



(16) 



P 



(A), nR{p,t)^ 



sin^ 7rQtop{t) 



P 



(17) 



Thus, for non-integer Qtop both A and B type fields produce an infinite number of gg-pairs 
with small momenta p 0. In the case of integer Qtop, Eq.(17) is not defined when p ^ 0. 
To obtain the correct expression for nR{p,t) in the above case, let us consider the system in 
a finite space volume V. Then, Eq.(13) must be changed as 



V/2 



nR{Pn,t) 



dxdye'P"^'-'-y^G^R\x - y) exp 



-v/2 



27r ^ ^R(kn,t) 



X ^Q-i>'n(t-x) _ ^ik„{t-x) _|_ ^ik„{t-y) _ ^-ik„{t-yY 



(18) 



Pm kn = 2Tm/V to be an integer. Furthermore, in the considered case the free Green function 
Go has the following form: 



G^S\x - y) ^ ^ j: e'^-^^-y^ 



Pn>0 



exp{2TTi{x — y) /V) 
V[l - eyi^{2'Ki{x - y) /V)] ' 



(19) 



Due to the {a+(p), a(p)} = V relation, to obtain the number of particles Nji{pm t) with the 
given momentum p„ we must divide nR(j)n,t) by volume V. Finally, we substitute instead 
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of in Eq.(18) its value at kn = equal to Qtopit) and obtain: 

V/2 

\ r QiPn{x-y)^2i{x-y)/V 



X exp 



Summing up the series in the exponential factor we find 

v/2 

Hence, for Qtop < we get: 

^VMPn,tj 1^ p„ = 0, 27r/V,...,27r(|Q,„,|-l)/r . ^^^^ 

We see that B-type field with an integer Qtop < produces exactly \Qtop\ of -R-quarks with 
vanishing momenta 0, 271/1-, 27r/V, — l)27r/V^. Simultaneously, the same field 

produces the equal number \Qtop\ of L antiquarks with identical momenta. For Qtop > 
quarks and antiquarks change places. Namely, in this case, Qtop of the L quarks and Qtop 
antiquarks are produced. 

The chiral properties of fermion systems in an external field were investigated in Ref.[9] 
just for the case of B-fields. The total chirality is changed according to the ABJ-anomaly 
[14] as 

K{t)-K{0) = 2Qtopit). (23) 

Here K{t) = J dxjl{t, x) is the total chirality at time t, j^{x) being the axial current density. 

The explicit calculations [9] for jQ{x,t) show that in the case of a finite volume system, 
there appear two different components in the density of chirality when a B-field acts. The first 
component stays finite at ^ — *■ oo. The second component contains terms ~ 1/V vanishing 
with V ^ oo but giving a nonzero result when they are integrated over V. These terms 
appear due to the discussed above phenomenon: particles with small momenta p ~ 1/V are 
necessarily produced by B-fields. So the corresponding part of the total chirality —2{Qtop) 
becomes delocalized over the volume V. The localized part being formed by contributions of 
quarks- antiquarks with finite momenta, is not conserved and has to obey just Eq.(23). But 
the total chirality (including the delocalized part) is conserved in the volume V!^ For V oo 
the described phenomenon represents the physical interpretation of the ABJ-anomaly for 
QED2 in the Coulomb gauge. 

^ It means that the equation for the divergence of axial current (ABJ-anomaly) changes for the case of 
finite volume as follows (Ref.[9]): ipnp{pn,t) + j{Pn,t) = — 3/27ri?(p„, i), n 7^ 0, j{pn 0,t) = K = 0. with 
{pn — 2TTn/V). So the total chirality being the particle number {Nr — Nr) — {Nl — Nj), is conserved for 
the fermion system in the external fields. 
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Moreover, this effect provides both the ABJ-anomaly and QED2 chiral condensates. The 
matter is the following. When is an integer number it is {Qtop\ of (IrQ_l{(1lQ.b) pairs to be 
delocalized with p ~ V^- ^^e case with an external field when Qtop is non-integer, in a 
coherent state with the distribution of Eq.(17) a non-integer number of quarks is delocalized. 
On the other hand, QED2 models have chiral vacuum condensates only when \Qtop\ is integer 
for essential self-consistent fields (see next Sec). This situation exists in the Schwinger 
model where \Qtop\ = 1 and the condensate is formed mainly by QRqLiQLQR) complexes. 

For the QED2 model with many electron flavours [4,6] there are no condensates without any 

( f) 

symmetry between flavours. But only in the case of symmetry the value Qlop becomes integer 
on important self-consistent fields for every fiavour /. The vacuum complexes exactly repeat 
the flavour symmetry in their structure (for example, a symmetrical condensate n^j^g^g{ or 

^f=iQRQL for the SU{N)-case). If there is no any flavour symmetry, there are no condensates. 
But chiral complexes with In V particles corresponding to Eq.(17) (/ dp/p ~ In V) are present 
here in vacimm states. We have quantitatively investigated these questions in Refs.[7,9]. 
Some properties of A and B flelds are very similar to instanton flelds in QCD [15]: 

i) Dirac determinant D{A) vanishes on the instanton held: 

ii) Euclidean winding number Qtop 7^ 0; 

iii) an instanton fleld also necessarily produces massless quark pairs [16]; 

iv) total chirahty is changed by the instanton fleld according to Eq.(23). 

But in the case of instantons D{A) vanishes due to the zero mode of the Dirac equation. 
Therefore, the Green function of a massless quark become s inflnite. That changes the quark 
chirality. The quark functions stay flnitc in 5-flelds of QEDg. For this reason, the chirality, 

in fact, appears to be the conserved number. EM- flelds can produce in QED2 only QhQr 
and ^L^L-pairs, but 2|(5top| of these particles appear immediately in the unobserved vacuum 
condensate. 

The essential flelds were studied for QED2 also in Ref.[12]. The consideration is based on 
the Euclidean functional integral method. Our consideration performed in the Minkowsky 
space and the Coulomb gauge to be used indicates an importance of EM fleld conflgurations 
different from those discussed in [17]. 

3 Physical role of the D(A) fields in Schwinger 
model 

We show in the present Section that EM- flelds where D{A) are responsible for 

A) both the absence of charged particles and for screening of charges; 

B) hadron states to have a deflnite chirahty {K — 0) in spite of the existence of the chiral 

condensate. 
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As it was stated in Sec. 2, Gauss integrals for matrix elements of the yS'(T)-operator are 
evaluated by a substitution of a saddle-point field A. This field depends on positions of 
particles in initial and final quark configurations. Those are the configurations determining 
investigated matrix elements. The formula for saddle-point field is obtained in Ref.[10] 
(Eq.(59) of Ref.[10]). One has to consider only small momenta p <^ m as the results of Sec.2 
witness. In the case of total charge Q = we have the following equation in the p — > limit: 



. g Ri{p) + Rf{p) cosm(T/2 - t) 
I ^ sinm-T/2 

.g Rf(p)-Ri{p) sinm(r/2-t) 
I ^ cosm-T/2 



(24) 



Ri{p) and Rf{p) are the sources depending on initial and final quark configurations: 



Rfip) = ^e(-p)(e-^^^'=-e-^K)+e(p)(e- 



(25) 



where x, x' {y, y') are coordinates of R{L) quarks and antiquarks in an initial state; x, x' {y, y') 
are those in a final state. The second term in Eq.(24) is inessential both for Green functions 
and D{A) because it vanishes being integrated over t. So it can be omitted. 

If charged particles exist, then charges at the large distance r (for example, q and q) 
inevitably would contribute to the wave function of two charged particles. Such S{T) matrix 
element will be independent on r within the r <^T time interval , i.e. within the region 

|xo - Xo| ~ l^o - XqI ~ T < r, |xo-Xo|~r. (26) 

But the moving {xq — > Xq) i?-quark contributes to the matrix element with the field strength: 



E^i\p, t) = — \Q{p) e-^P^° - e(-p) e-'f^° 

TT}j 

or in the coordinate space: 

E^'^Xx^t) 



2nm 



1 



1 



X 



Xq + iO X — Xq — tO 



cosm(T/2 — t) 
sinm • T/2 



cosm(T/2 — t) 
sinm ■ T/2 



(27) 



(28) 



This is just the field of the A-typc. To obtain the field strength for an i?-antiquark we take 
— >■ Xq, Xq — > Xq and change the general sign in Eqs. (27) and (28). 

Substituting Eq.(27) into Eqs. (7), (9) we obtain that the gnqR configuration, contributing 
to S{T) as (Eqs. (3), (9)): 



bR{x'Q)aR{xQ)G^P (T, 4; 0, x'q) D{A(^+9^)G^P (T, xq; 0, x[,)a+(xo)6^ (x[,) 



(29) 
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decreases at r — > oo as r ^ (^D(A^'^^) ~ r^, G ^ r ^j. Let us note that D(A'^'^^) does not 

vanish but increases with r ^ oo. This is a direct consequence that E^i^ is a purely imaginary 
value at {xq — Xq) of Eq.(28). However it is shown in the Appendix that the result can 
be obtained also by integrating over the following real fields of the A-type: 



E{p,t) = i(3 Q{p)e-'P^° - e(-p)e" 



■ipxo 



cosm{T/2-t) 
smm-T/2 



(30) 



where /9 is a real parameter. The r dependence in the discussed amplitude is due only the 
integration over the /5 parameter in (30). In this case wc have that D{A) ~ l/r^'^ and 



Gr,l ~ the Jacobian being 

gives the same result as Eq.(28): 



'Inr. Therefore, the integration over fields, Eq.(30), 



1 1 



at 



oo. 



(31) 



Thus, the screening mechanism is based on the fields of the class A with D{A) — > 0. 
It means that S{T) contains only configurations with q and q near each other. Additional 
particles should be present in configurations if some quark situated far from some antiquark. 
Their role is to compensate charges of distant qq particles. We have seen in the previous 
section that the smallness of D{A) can be compensated by a contribution of quarks with 
small momenta p ~ 1/r. One can demonstrate this phenomena for the considered case. For 
this aim let us substitute the EM field, Eq.(27), into the factor 



/ 



dx dx'G^^\x, T, x', T)aJ(x)6j(x') 



(32) 



representing in S{T) an additional qjtqR-peiiY in the final state. If the pair creation happens 
near the quark xq 



Xq, we can neglect the infiuence of the antiquark (xg — Xq ) field: 
X ^ \x — Xq\ — \x' — Xq\ ~ |,Xo — Xq\ ~ Tq. Using (7), we obtain the Green function 
G'j^{x, T; x', T) for the field (27). In the region |a; — Xq] 3> |a;' — Xq] ~ T this Green function 
becomes: 



GS,^)(T,x;T,x') = -^exp 



1 



2m 



± 

I 



dti 



cos m(T. 2 — ti 
sinm • T/2 



In- 



1 



T — ti — X' + Xq — iO 



(33) 



In evaluating Eq.(33) we use the following formula 

„(o) / „ \ 1 1 A ^ \ f dxidti 
G)i'{x,T;xi,ti)dxidtiAQ[xi,ti) = j 

dxidti 



Ao{xi,ti 



2Tii 



2m T — X -{- xi — ti — 
ln(r-x + xi -ti)E{xi,ti) ■ 



(34) 



The Green function (33) does not depend on the quark coordinate x. It depends only 
on the antiquark coordinate x' . Thus the additional gi^g^-pair in the wave function (29) 
represents a quark with the zero momentum (/ a%{x)dx — a^{p — Q)) and an antiquark with 
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a finite momentum determined by the field strength E^^\x,t). The quark x is situated for 
from the original one {xq ~ Xq), while the antiquark x' is near xq, and, hence, its charge 
participates in the screening. The same effect happens with L-Green functions. 

Thus A-fields prohibit the existence of charged configurations and inevitably provide the 
last ones with local screening charges. 

Now we proceed with the B-type fields. They maintain the chirality conservation of any 
quark configurations. In order to prove the statement let us consider an appearance of an 
additional pair a^(a;o)&^(yo) the final state, i.e. the changing of chirality by two. Eq.(24) 
shows that such a pair produces the B-type field: 

t) = \e{-p)e-'P^° + 0(p)e'^™^l ^ (35) 



m 



sinm(T/2) 



QSpiT) = -|- £ dt E(^)(0, t) = l. (35') 



The B- field nullifies D{A) and prohibits the suggested process. It generates additional 
gg-pairs. Indeed, let us consider the factor representing two such pairs: (iR{x)b~^{x') and 
(it{y)^L{y')- The wave function acquires the factor 

I dxdx'dydy'aUx)GPix,T;x\T)bUx>tiy)GP{y,T-,y',T)bt{y') • (36) 

For the fields (35), the Green functions Gr and Gl do not depend on x and on y' respectively 
Thus E^^\p^t) produces a chiral pair qRqi near the original chiral pair and g/^g^-particles 
with small momentum to be far from the region E(^) 7^ 0. The first chiral pair compensates 
the change of the local chirality carried by an the original one. The existence of the chiral 
condensate seems to be very natural when the creation of chiral pairs with vanishing total 
momenta appear. As it was explained in Sec. 2, condensates really arises in QED2 model when 
Qtop is an integer number for all essential B-fields. Such conditions exist in the Schwinger 
model (see Eq.(35')), and in the many electron model [4,6] with some symmetry between 
flavours. 

Let us note that B-type fields prohibit the change of chirality, but states with the definite 
chirality can exist. For example, the configuration with a QrQl pair both in the initial and 
in the final states produces the field: 



Eix t) = cosm(r/2-0 
^ ' ^ 27rm sinm-T/2 



+ 



X — Xq + iO X — y'd — iO x — Xq — iO x — yQ + iO 

(37) 

For this field D{A) is non-zero and additional pairs need not be produced. These quark 
configurations form the chiral condensate in the Schwinger model [7]. 

To summarize, we see that localized states only with the definite chirality can be station- 
ary in QED2. The chirality exchange between them and the vacuum condensate is impossi- 
ble. As we have already noted that is the property which makes possible the confinement of 
fermion numbers in QED2. 
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4 Space-time picture of charge screening and of the 
hadronization in physical processes 

The production of a chiral quark pair QRqLiQR, Ql) by an external source is an analogy of the 
e'^e~ hadron process in QED2. The process had been studied in Refs. [5,8,9]. The formulae 
for EM charge densities for R- and L-particles were obtained there to be 



2 ^--x 



(38) 



Ji being the Bessel function. One can observe in Eq.(38) the leading and the screening 
charges. The both arc situated near the light cone. The screening charge is accumulated in 
the region Ax ~ 1/rrPt near the leading particle. After the time t ~ pin/m^ {Ax ~ and 
X ~ in the restframe) the leading quark starts to interact with the screening ones. The 
leading quark moves as a free particle till times t < Pin/rn^ [5,8]. 

We consider in this section quantitatively all three stages for evolution in time of the 
parton process: a) neutralization of jets, b) accumulation of partons, c) transition into 
hadrons. The physics is closely connected with the field generated by a qnqL pair, i.e. with 
the B-type field. 

Let us formulate the results. 

1. R and L jets practically do not effect with each other after the time to ~ m~^, when 
particles to be already generated neutralize total charges and chiralities of each jet. 

2. The neutralizing charges appear as a chiral pair [quqi) with momenta p m. This 
effect is the result of B-field action. 

3. Quarks of large momenta p ^ m are accumulated in almost neutral systems {t > to). 
They exist only at the hght cone region. They softly turn into hadrons, which slow down 
~ p, the hadronization time is i ~ p/m^. 

4. Inside the light cone there are only quarks to have been organized into hadrons. In 
the region between charges there are no any string of quarks- ant iquarks. 

The screening of total jet charges {qji and q^) and the transition into neutral hadrons 
require to change chirality of the g'^g^-pair {K — 2). As we know (Sec.3) the field of B- 
type must appear with such effect. Indeed, calculating the field strength E{x, t) for the 
distribution of Eq.(38) we obtain 

E{x, t)=gJo {Vf" - x^) e{t^ - x^) . (39) 

It is the field of the B-type since at p — > 

^, , ^ sin mi 

E{p,t) = 2g . (40) 

m 

The corresponding shift of chirahty is equal to 

K{t) - K{0) = 2Qtop{t) = 2{cosmt - 1) . (41) 
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The chirality changing (41) imphes the creation of two kinds of quark pairs by the B-field 
(see Sec. 2-3). 

1) One additional quqi pair with the momentum p m. These particles compensate 
total charges and chiralities of R and jets. This pair exists after to ~ 1/m. 

2) One gnQL pair with the total momentum p ~ 1/i. It is expanded always over the 
whole causal region {2t). This pair appears simultaneously with the first one, as chirality 
of local QED2 systems can be changed only on the time of fluctuations (~ 1/m). Chirality 
oscillations mean a fluctuative exchange between just this pair and the condensate. The pair 
joins the condensate at t — > 00. 

The latter pair has no connection with the parton accumulation and transitions into 
hadrons. The both processes generate independent and neutral {Q — 0, K — 0) compact 
R, L packets. The quarks of the first gnqi pair participate in the hadron production. They 
form hadrons with p ~ m at the time t 1/m. 

To prove the above assertions, we investigate the time dependence of the number of 
quarks produced by the field, Eq.(39), in the considered process. The strict calculations by 
means of the bosonization method gives [9] the answer very similar to Eq.(13). The density 
of the right-handed quarks is 



Here, a{x,t) can be expressed in terms of the field strength E{x,t), Eq.(39), by means of 
Eq.(14). The function 



determines the mean density of the quarks in the vacuum. 

The integrand in Eq.(42) depends on ^ = t—x and ^' = t—y. For 1/m andp < m 

only the difference ^ — ^' is restricted by a rapid decrease of the function f{x — y) — /(O- 
The values of ^ and ^' themselves are large {—t). In this case the integration in Eq.(14) goes 
over all possible x\ and ti. Thus, we have 



f dxdy e^v^^-y) 



(42) 





(43) 



[x^t) ^ -^Q{t - x) [ dti I gJo{m){tj- xj)Q{tl- xl)dxi 



Jo 



(44) 



Applying Eq.(16) and introducing the variables C — C ~ C and we obtain 




(45) 



where 



F(C) = Cexp[z7iQtop{t)e{C)]2tsm2nQtop{t) , 



(46) 
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and Qtop{t) is the same as in Eq.(41). Finally, substituting Eq.(46) in Eq.(45) we get 



nnip, t) = A{p) sin^ 7rQtop{t) + B{p) sin 27rQtop{t) , (47) 

and 

Mp) = - / ^C/(C) sinpC, Bip) = - / dCm cospC ■ (48) 

TT JO TT Jo 

Equation (47) clearly demonstrates the oscillating nature of the density of quarks produced 
with small momenta. The first term in Eq.(47) describes production of quarks by the field 
given by Eq.(39). For /(C) = e~'^^ {5 — 0) wc return to Eq.(17). The second term describes 
the influence of the field on vacuum quarks. That is the effect which produces oscillations. 

R quarks of Eqs.(42)-(48) are bounded into chiral pairs. Since in Eq.(41) there are no 
restrictions on x and y values (|a;|, \y\) < t), the pairs are expanded over the whole volume 
~ 2t, i.e. the total momentum of the pair is p ~ t~^. The momenta of the quarks inside one 
pair are determined by the integrals in Eq.(48) and are of the order of p ~ m. The total 
number of bounded quarks is 

/ ^ / ^^^^^ ^ ^^^^ 

Therefore we have dealt with one pair. The value of Qtop oscillations also confirms this fact. 

Let us now consider particles with momenta p ^ m. They do not participate in transi- 
tions into vacuum, but just these particles form hadrons. The last ones with momentum P 
are formed during the time 

i ~ ^ . (50) 

The size of the quark packet, Eq.(38), becomes then equal to the size of such hadrons ~ 1/P. 
After the time ~ P/m^ all quarks with momenta p < P have to be bounded. Equation (42) 
gives the opportunity to verify this qualitative picture. In some aspects we repeat Ref.[9]. 

For particles withp ^ m the value of |a; — 1/| becomes ~ 1/p, and f{() ~ /(O) = 1. These 
quarks do not interfere with vacuum quarks. We apply again Eq.(16) and rewrite Eq.(42) 
in the form {( — {t — x)): 



Jp ZTT J-oo 



2 



(51) 



Since jets become neutral at t ~ 1/m we have = 1%^ = ni = n^. 
In Eq.(14) the essential integration region is near to the light cone: 

0<(i^ti-Xi<C^t-xr^-. (52) 

P 

Therefore, quarks with momenta p are produced only in the region 1/p near the leading 
particles. Inside the light cone there are quarks with p < mH. They have already formed 
hadrons. 
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According to Eq.(52) the — x\ interval is approximately t{ — x{ — 2tiC,i, and a is 

.2 . 



«(C, ^) = / '^^i I dCiJ (m^^ = e(C) [jo (m^C) - 1 



(53) 



since Q{t\ — x\)Q{C, — ti + xi) = 0(Ci)6(C ~ Ci)- Eq-(53) shows that jets are independent 
from each other for {t > 1/m) because only the field of the i?-jet produce i?-quarks. In the 
time interval <^t<^ p/m^ we can expand the Bessel function in Eq.(53) and get 

nR{p,t)^-[ -, m-^<t<-^. 54 

& \ p ) p 

At times t > p/m^ the argument of the Bessel function becomes large and we expand the 
exponential factor in Eq.(51) in powers of a{C,t). The integrations over k and ( give the 
distribution to be independent of time as follows 

nRip,t)^—, t>^. (55) 
p 

At times t ~ p/rn^ the distribution of Eq.(54) becomes the order of unity. Quarks with 
the momentum p begin to turn into hadrons. Eq.(55) represents the distribution of quarks 
inside hadrons. 

Indeed, we have in the parton model 

r dP 
nnip, t)=j N{P, t)nR{P, p) — . (56) 

Here N[P, t) is the number of hadrons with the momentum P. In the process considered 
this quantity is equal to [5,9]: 

^ ' <*(«) I *W) ^ ' 

Here, A^{p) are hadron operators. The distribution of quarks inside a hadron with the 
momentum P is equal to: 

nR{P,p) = (*,(P)|a+(p)aK(p)|*,(p)) = ^ . (58) 

We obtain Eq.(58) by means of the parton wave function for a QED2 hadron. We had found 
in Refs. [9,10] that 

/27r f dk ^ 

m 



/ zTT c nk' 

^^HiP) = y ^ / ^a+(P - k)h^R{k)\n, > , P » 

aR,L(A;)|Jlo >= &r,l(A;)|^^o >= , for A;> m . (59) 

Substituting Eqs. (57)-(58) in Eq.(56) we obtain the distribution of Eq.(55). 

The leading particles lose their energies and turn into hadrons after the time t ~ Pin/rri^. 
Hence, the most rapid hadrons are produced in the last turn (it is similar to Refs. [5, 17]). 
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5 Conclusion 



In the present paper it was demonstrated that confinement of the charged states in QED2 
is provided by the field configurations which nuUify the Dirac determinant while Green 
functions stay finite in these fields. It is possible that this confinement mechanism can 
exist not only in two dimensions. But at present it is unknown whether the fields with 
similar properties play any role in QCD, in four dimensions. But if the D{A) ~ fields 
are important also in QCD, the mechanism based on such fields would be very attractive 
because it naturally provides theory with many phenomenological inviting properties. 

In fact, the colour screening in QCD could appear because quark-antiquark create a 
coloured field providing for the D{A) determinant to go to zero when the distance between 
them increases. This field produces (with probability of one) soft quarks which screen the 
initial colour. These quarks create the baryon number current which could, generally, screen 
also the triality of the sources. 

Certainly, the screening is much more complicated for QCD because of its non-Abelian 
nature (also we must screen the gluon colour) and probably needs more time than in QED2. 
However in the described scenario the momenta of leading quarks are much more than 
momenta of the screening ones. So it seems that the time of screening would be less than 
the time of hadronization. As a consequence hadrons are produced after the time when the 
screening of colour is over and jets are independent and colourless. Thus the fast hadrons 
(which are produced the last) would be created quite independently in different jets. Light 
quarks would play a fundamental role in this confinement mechanism. 



The explanation concerning the imaginary saddle field in Eqs.(28),(29) consists in the explicit 
evaluation of the functional integral. We have for the matrix element considered: 

So = SlliT) = j DA{x) eMiS{A)]G^P (T, 4; 0^) D(A)g''P{T, xo; 0, xo) . (60) 

We are interested only in components of A with very small p, these ones are important for 
the dependence on the distance r between q and q when r — > 00. Then we have 

SiA) = / ^ |!!!^a(p)a(_p) + ^g[R,(p) + ^R,{p)]a{-p)\ + S{p » po) , (61) 

where ^ 

Po<m, a{p)= [ Ao{p,t)dt . {A.2') 
Jo 

We neglect in Eq.(A.2) the kinetic energy of the small components: Jg^" p^Aq{p, t)dp/27T. We 
take the sources from Eq.(25) and rewrite down them for the qq-case in the following way 

Riip) + Rf{p) = Riip) + R2{p) . (62) 



16 



Ri{p) — Rl{—p) and R2ip) — —R^i—p) are real and imaginary parts of the sources. 
The integral can be easily transformed to 

So^ZxJ d(3e''^^ I Da{r)5 ( J Ri{p)a{p)^ - (3x\ 



exp < I 



dp 
2^ 



im'^\p\ 



a{p)a{-p) + igR2{p)a{-p) 



GrGrD 



(63) 



Here Z includes all contributions independent on r, and 

^-1 



Riip)Rii-p) dp 

~ m r . 

2n 



\P\ 



Let us change the integration variables according to 

a{p) = + d{p) = ao{P,p) + dip) . 



(64) 



P 



We obtain 



So^ZxJ dpD{ao)GR{ao)GR{ao) ■ I 



(65) 



(66) 



/ = /Da(p)5(/i?.a|)exp[-/ 



dpm'^\p\_, 
^^^a{p)a{-p) 



The dependence on i?2 is inessential for r oo. Thus the integration in Eq.(47) goes over 
real fields and these fields are of the A class. Practically we finish the proof of our assertion 
in the text. The further evaluation of / is straightforward. 



I = 



/D5/|exp -/ 



dp TTi^ \p\ 

— ^^d{p)d{-p) + iCRiip)d{-p) 



dC 
2ti 



exp 



^Inr 



(67) 



Substituting all factors in Eq.(A.7) we get the result of the text. Obviously, the saddle point 
for integration over (3 is imaginary: (3 = i, and we return to the field of Eq.(22). 
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